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Using the CLEO II detector, we have measured the differential cross sections for exclusive 
two-photon production of light pseudoscalar mesons 7r°, 77, and 77'. From our measurements 
we have obtained the form factors associated with the electromagnetic transitions 7*7 —> 
meson. We have measured these form factors in the momentum transfer ranges from 1.5 to 
9, 20, and 30 GeV^ for tt", rj, and 77', respectively, and have made comparisons to various 
theoretical predictions. 
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I. INTRODUCTION 



Production of even C-parity hadronic matter in e^e~ scattering provides a unique opportunity to study the 
properties of strong interactions. To leading order in quantum electrodynamics (QED) these processes are described 
as the interaction between two photons emitted by the scattered electrons^ Although in e^e~ scattering the probe 
and the target are both represented by photons that are carriers of the electromagnetic force, these space-like photons 
can produce a pair of quarks that interact strongly and are observed in the form of hadrons. Therefore, by measuring 
the four-momenta of the scattered electrons we can study the dynamics of strong interactions. The quantities of 
interest in these studies are the form factors associated with the transitions between the photons and the hadrons. 

This Article describes the measurements |l| of the differential cross sections for the production of a single pseu- 
doscalar meson in e^e~ scattering: 

e+e" e+e"7^, (1) 

where 7?. is a tt", r; or 77'. We measure these cross sections in a "single-tagged" experimental mode where one 
of the scattered electrons is detected ("tagged"), while the other electron is scattered at a very small angle and 
therefore remains undetected ("untagged"). The mesons produced in e^e~ scattering are observed through their 
decays to various fully reconstructed final states. The tagged electron emits a highly off-shell photon (7*), whereas 
the untagged electron emits a nearly on-shell photon (7) . We measure the dependence of the meson production rate 
on the squared momentum transfer carried by the highly off-shell photon. This momentum transfer is determined 
by energy-momentum conservation as applied to the tag: 

= -{pi,-ptf ^2EbEt{l- cos Ot), (2) 

where pt, and pt are the four-momenta of the incident beam-energy electron and the tag, and Et are corresponding 
energies, and 9t is the scattering angle^. From the measurements of the differential rates 

da{e^e~ — > e~TZ) , , 

we obtain the transition form factors J-~f'~,Tz that describe the effect of the strong interaction in the 7*7 —* TZ 
transition amplitudes. 

To relate the differential cross sections to the transition form factors we employ the theoretical framework developed 
by V.M. Budnev et al. (BGMS formalism). In BGMS the process e^e~ e^e~TZ is divided into two parts: 
e^e~ — > e"^e~7*7 and 7*7 TZ. The first part is completely calculable in QED and the second part is defined in 
terms of the transition form factors ^niO^) ■ In the case of pseudoscalar mesons there is only one form factor. At 
zero momentum transfer this form factor is expressed as: 

\T m^P- ^ 64^^(7^^77) 

where a is the QED coupling constant, Mu is the mass and T{TZ 77) is the two-photon partial width of the meson 
TZ. The transition form factors cannot be calculated directly from Quantum Chromodynamics (QCD). However, they 
have been estimated using perturbative QCD (PQCD), a sum-rules approach, and other theoretical methods. 

One of the important concepts of PQCD-based methods is a factorization procedure that separates perturbative 
short-distance effects from non-perturbative long-distance ones. While the former are understood well and can be 
calculated using PQCD, the latter are known only asymptotically, in the limit 00. In PQCD-based calculations 

the transition form factor T^^-^-n is expressed as a convolution of a perturbative Hard Scattering Amplitude (HSA) |^ 
and the soft non-perturbative wave function of the meson. 

Brodsky and Lepage employed PQCD to find the asymptotic behavior of the 7*7 TZ transition form factors in 
the limit ^ 00 |||: 

lim QV-y-^ii(Q') = 2/7J, (5) 

where f-n is the meson decay constant. In addition, it has been predicted that in this limit any mesonic wave function 
evolves to the asymptotic wave function of unique shape . 



^Unless otherwise specified, we use the term "electron" for either an electron or a positron. 
^ The electron mass is neglected in Eqn. ^ 
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While PQCD predicts the form factors of the 7*7 ^ TZ transitions at large momentum transfer, the behavior of 
these form factors in the limit ^ can be determined from the axial anomaly in the chiral limit of QCD. 
For -K^ and f] the axial anomaly yields 

lim Tj'^niQ^) = . Ir , (6) 

to leading order in m^/M^ and m'^/M'^ where rriu and nid are the masses of the u and d quarks. This prediction 
does not hold with the same precision for rj' due to the larger value of the s-quark mass. In addition, even if the 
s-quark mass were small, this prediction might be broken for rj' because this particle is an unlikely candidate for the 
Goldstone boson [^,0. 

To describe the soft non-perturbative region of a simple interpolation between — > and — > 00 limits has 
been proposed 

To quantify the long-distance effects in the soft non-perturbative region, Chernyak and Zhitnitsky employed the 
sum-rules method IllJ to derive the wave function of the pion at experimentally accessible momentum transfers (the 
CZ wave function) W. They demonstrated that the proposed wave function successfully describes experimental data 
on the Xc decay into two pions and the electromagnetic form factor of the charged pion. However, because the 
theoretical predictions for these processes depend on the strong interaction coupling constant Qs, this introduced a 
large uncertainty in the determination of the CZ wave function. 

Since the asymptotic and CZ wave functions were proposed, they have often been used to describe the non- 
perturbative parts of transition amplitudes in various PQCD calculations. Jakob, KroU, and Raulfs employed these 
wave functions and PQCD to calculate T^.^^o |^,^. These authors have also taken into account small QCD 
radiative corrections, incorporated into the PQCD technique by Lee and Sterman KroU has concluded that the 
CZ wave function disagrees with our preliminary results On the contrary, a competing perturbative analysis 

of Cao, Huang, and Ma jl^ yielded that either the asymptotic or the CZ wave function is sufficient to describe the 
data. These authors took into account quark transverse momentum corrections and neglected the QCD radiative 
corrections, estimating the latter as small. 

While PQCD-based methods are often employed to predict rates for exclusive processes^, the applicability of these 
methods at experimentally accessible momentum transfers remains one of the outstanding problems of the theory of 
strong interactions. Extensive discussion of the validity of the PQCD approach can be found in the literature Ipij-pst. 

To avoid ambiguities of the PQCD-inspired calculations at of the order of several GeV^ , Radyushkin et al. devel- 
oped an approach ]2e|"|3i| based on the sum-rules method [|ll| that they employed to predict the 7*7 — > tt" transition 
form factor This prediction depends on the model of the hadronic spectrum chosen to describe an almost real 
photon emitted by the untagged electron. It also depends on the values of vacuum condensates which represent 
non-perturbative matrix elements. The theoretical result of Radyushkin et al. reproduces the PQCD-predicted 
shape of the transition form factor but disagrees with the absolute value given by Eqn. m by about 15% in the limit 
00. The authors have stressed that this discrepancy is irrelevant in the region of Q below 10 GeV^ and could, 
in principle, be eliminated by including the QCD evolution into the theoretical analysis . It should be noted that 
the discussed theoretical analysis exactly reproduces the asymptotic prediction of PQCD given by Eqn. ^ when both 
photons are highly off-mass shell. We should emphasize that at present the non-perturbative treatment of various 
exclusive processes in a way similar to the approach of Radyushkin et al. is the subject of significant theoretical 
interest. For example, the QCD sum-rules method has been employed recently to predict the form factors in the 
semileptonic decays of the B mesons^ p4|-|3^. 

The 7*7 TZ transition form factors have been studied by several experiments. The LEPTON-G experiment 
measured T-^*^n and T.^*^ri' in the time-like momentum transfer region up to 0.24 GeV^ using the rare electromagnetic 
decays -q — > fi'^fj,~'y and rj' — > ji^ ji~'y [^. In order to achieve higher values of , the space- like photons produced 
in two-photon interactions were utilized by the PLUTO experiment to measure Tr^i-^i^i up to 1 GeV^ @ and by 
the TPC/27 collaboration to study J-~i»~ir, and J?%y..y^' up to 7 GeV^ [^. More recently, the CELLO experiment 
measured T.^t..^.^o at up to 2.7 GeV^ and T.^-'-.^n and T.^-„i^i at up to 3.4 GeV^ jsof . 



For example, these methods have been utilized to calculate the nucleon form factors [0jl8| and the 5° 
branching fraction see also [^o|J^. 

* Recent results of other theoretical developments relevant to our experimental study can be found in the litera- 
ture |||-||. 
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We employ two-photon interactions to measure the transition form factors J-~,*~fn in the space-hke regions of the 
momentum transfer between 1.5 and 9 GeV^ for -k^ , 1.5 and 20 GeV^ for 77, and 1.5 and 30 GeV^ for r}' . We study 
the transition form factors of tv^ , rj, and 77' using the decays: 

71-0 _^ 



V 77, 
r; ttOttV" ^ 67, 

T]' —> pO-y _^ 7r^7r^7, 

' + — + -o 

rj n TV rj ^ TV TV Z7, 

T]' —^ n'^Tv'^ri 67, 
T] -K^-K^ri 57r" IO7, 



T]' ttOtt"?; Sty^'tt^tv — > tt^tv 67, 
rj' Tv'^'K~ri —* 7r'''7r~37r" ^ 7r^7r~67. 



We analyze the last two decay chains of rj' together since they are observed in the same final state 7r'*'7r~67. 

This paper is structured as follows: Section II describes the CLEO II detector and the data sample that we use 
for our measurements. Event selection criteria, experimental technique, and the analysis procedure for 77 final 
states are explained in Section III. Analyses of other final states with only photons are described in Section IV and 
analyses of final states with charged pions are described in Section V. The unfolding procedure for the transition form 
factors is described in Section VI. The results are compared with some existing theoretical predictions in Section VII. 
Conclusions are presented in Section VIII. 



II. EXPERIMENTAL APPARATUS AND MONTE CARLO SIMULATION 



A. The CLEO II Detector and Data Sample 



The CLEO II detector |5l| is a general-purpose magnetic spectrometer which provides energy and momentum 
measurements for elementary particles. It is operated at the Cornell Electron Storage Ring (CESR), a symmetric 
e^e~ collider running at a center-of-mass energy near 10.6 GeV. The major objectives of the CLEO experiment are 
the studies of the properties of heavy mesons that contain b or c quarks. However, owing to the versatility of the 
detector, analyses of tagged and untagged two-photon interactions, detailed studies of r-lepton decays, and careful 
examination of quark and gluon fragmentation and other processes are also possible. 

The active components of CLEO II include central tracking detectors, time-of-flight (TP) scintillator counters, 
muon detectors, and a Csl calorimeter for electromagnetic showers. The calorimeter consists of a barrel part covering 
polar angles above 37° and two endcap parts each covering the region between 13° and 37°, where the polar angle 
is measured with respect to the beam axis. The energy resolution of the barrel calorimeter for photons of energies 
above 500 MeV is 2%. The central tracking detectors consist of three concentric cylindrical drift chambers that cover 
the polar angles above 18°. From smallest to largest radii these are: the precision tracking layers detector, the vertex 
detector (VD), and the main drift chamber. The measurements of the specific ionization energy losses in the outer 
layers of the main drift chamber and flight times in the TF system provide discrimination between charged particles 
of different species. All detector subsystems except the muon detectors reside in a uniform axial magnetic field of 1.5 
Tesla. 

The data sample employed in our analysis corresponds to an integrated e^e^ luminosity of 2.88 ± 0.03 fb^^. Two 
thirds of the data was collected at e''"e~ center-of-mass energy of y/s — 2Eb = 10.58 GeV, the remainder at 10.52 
GeV. 
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B. Trigger System 



The CLEO II detector has a three-level hardware trigger system followed by a software filter. The fastest, 
"zeroth" level (LO) trigger can be either track-based (using the VD and TF) or energy-based (demanding a minimum 
energy deposition of about 500 MeV in the Csl calorimeter). The calorimeter LO information develops slowly, so it 
only comes into effect if the track-based LO trigger fails; in such cases the tracking information is lost. 

The first level (LI) trigger uses track-based information from the VD, TF, and main drift chamber; tracks of 
transverse momenta in excess of about 340 MeV/c are identified by either of two independent track processors 
employed in the trigger decisions. Calorimeter information is also utilized at LI. High threshold bits, designed to be 
set by showering particles, have a threshold of about 500 MeV; low threshold bits, designed to trigger on minimum 
ionizing particles, have a threshold of about 100 MeV. To trigger at LI on two low-energy clusters, they must be 
well-separated in space. 

More detailed information from the VD and main drift chamber is used in the second level (L2) trigger. The 
requirements and accessed momentum range varied between data subsets, but are all modeled in our detector simula- 
tions. The software filter (LVL3) is optimized to suppress backgrounds from interactions of the beams with residual 
gas and vacuum chamber walls. Events which pass LVL3 are recorded. In addition, every eighth event that fails 
LVL3 is also recorded to allow the LVL3 efficiency to be studied. 

The efficiencies of the various trigger components have been measured using data collected with independent or 
partially independent simultaneous trigger requirements and are incorporated in the detector simulations [Q,^j5^]. 
The simulations are carefully run to match the integrated luminosity associated with each trigger configuration. This 
is necessary because exact trigger requirements in CLEO II have been changed over time to improve the trigger 
efficiency for events of low particle multiplicities. 



C. Monte Carlo Simulation 

In our analysis we use a two-photon Monte Carlo (MC) simulation program j55| that is based on the BGMS 
formalism H. The 7*7* TZ transition form factors are approximated by: 



\P'y''y*Tz{Q ,q)\ = \J^r'-in{Q 



647^^(7^ 77) 



(47ra)2 Ml {1 + / h.\Y {I + / h.\Y 



(8) 



where and are the absolute values of the squared four-momenta carried by the space-like photons. The pole-mass 
parameter Km = 770 MeV has been chosen to approximate the momentum transfer dependence of the form factors. 
It should be noted that while we have chosen this parameter to be practically the mass, as predicted by the Vector 
Meson Dominance (VMD) model the pole-mass behavior of the transition form factors J-~f»^u and the value of 
the parameter K-ji in the range between 700 and 900 MeV are indicated by various theoretical predictions p^ , |57{ 
that are not based on VMD. Notice that in the approximation given by Eqn. |^ we assume a factorization of the form 
factor into the q^- and Q^-dependent parts |^^. In the same two-photon simulation program we also generate the 
decays of the produced mesons. To account for the relativistic effects, helicity conservation, and presence of spin-one 
particles we simulate the decay chain r)' p^-y — > 7r^7r~7 according to: 

d^r(V^p°7^^+^-7) ^ ^.^2 m,r{m^^) 

d cos 6* drn^rr m-^-n {raj, — mi-^Y + rajtV^ira-n-K)^ 

with the energy-dependent width, r(m.^„), parameterized by: 

r(m,.) =r(mp)^, (10) 

\Pn\ 

where 6* is the angle between the directions of one of the charged pions and the signal photon, i?-,, is the energy of 
the photon, T{mp) — 151 MeV and mp = 768 MeV/c^ are the nominal width and mass of p" |59| , m^^ is the actual 
mass of p", and and |p7i| are the magnitudes of the charged-pion momenta for the actual and nominal masses of 
p°, respectively. The charged-pion momenta and the angle 6* are defined in the center-of-mass frame of p" and the 
energy of the photon is defined in the center-of-mass frame of 77'. 

The transport of the generated MC particles through the CLEO II detector is performed by a GEANT-based 
detector simulation program. The generated events are then processed by the event reconstruction program which also 
"simulates" random electronic noise and beam-related spurious energy clusters by adding hits from random-trigger 
data samples into the MC events. 
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III. ANALYSES OF SINGLE-TAGGED 77 FINAL STATES 



A. Trigger 

The single-tagged two-plioton reactions e+e" e'^e~n° and e+e" e^e'rj followed by the decays 7r°^77 and 
r?— >77 are recorded using either a track-based or an energy-based LO trigger. The track-based LO trigger is satisfied 
when the scattered electron passes through the VD and enters the endcap TF. For tags that scatter at polar angles 
above 24.5°, thus passing through the entire VD volume, the efficiency of this trigger is about 80% and is determined 
by the size of the wire-chamber drift colls compared to the time allowed to make the LO decision. At smaller polar 
angles we rely on the energy-based LO trigger. The efficiency of this trigger is 98% (100%) for electrons which 
deposit 1.0 GeV (more than 1.6 GeV) of energy in the calorimeter. 

The LI trigger is satisfied when at least two clusters, each of energy above 500 MeV, are detected in the calorimeter, 
with one in the barrel region and the other in one of the endcap regions. There are no L2 requirements for events 
passing the LO and LI trigger conditions described above. 

To be recorded, events must fulfill the transverse-momentum requirement of the LVL3 filter that assigns momenta 
to all calorimeter clusters assuming that they are photons produced at the primary interaction point in the geometrical 
center of the CLEO II detector. This LVL3 criterion rejects events if the net vector momentum has a component 
normal to the beam axis in excess of 0.7 GeV/c (1.4 GeV/c) when the total energy detected in the calorimeter is 
larger than 1.0 (5.0) GeV. 

B. Analysis Procedure 

In the first part of this section we describe the event selection criteria based on the event topology for the signal 
production processes. In the second part we explain selection criteria aimed at the suppression of random background. 
In the third part we discuss the event quality requirements designed to isolate signal events with large uncertainty 
in the detection efHciency. Finally, in the last part of this section we show the invariant mass spectra for data events 
that fulfill all selection criteria. 



1. Basic Selection Criteria 

The event selection criteria for single-tagged 77 final states are designed to isolate two-photon events for which 
the trigger efficiency is high and in which the only missing particle is the untagged electron of high momentum. 
These events are characterized by the high-energy shower produced by the tag in the endcap calorimeter and two 
electromagnetic showers of total energy larger than 1 GeV produced by the photons in the barrel calorimeter. 

We select events in which three or four energy clusters and no more than one charged track have been reconstructed. 
The energy of each barrel (endcap) cluster must be larger than 30 (50) MeV. The most energetic cluster is assumed to 
be produced by the tag and must be in the endcap calorimeter. If a charged track is found, its projected intersection 
point with the calorimeter must agree with the tag's shower position within 20° as estimated at the primary interaction 
point. The position of each shower is determined from the energy- weighted average of the centers of the crystals 
forming this shower. To provide an efficient trigger, the energy of the tag candidate detected in the calorimeter should 
be above 1.0 GeV (at a later stage of the analysis procedure this cut will be superseded by a tighter requirement). 
Out of the remaining energy clusters, the two most energetic must be found in the barrel calorimeter at polar angles 
above 45° {i.e., excluding calorimeter edges), and are assumed to have come from the -k^ or r) decays. The fourth 
energy cluster, if found, should contain less than 200 MeV of energy; the efficiency loss due to this requirement is less 
than 0.25%. Events with this additional energy cluster may be either signal or beam-gas events with a beam-related 
noise cluster or partially reconstructed background events of higher particle multiplicities that mimic single-tagged 
tt" or 77 production. By allowing an extra energy cluster, we reduce the uncertainty in the signal efficiency while 
providing the opportunity for background estimates. A tighter cut on the energy of an additional cluster would make 
our results more sensitive to the modelling of the noise-related energy clusters and a looser cut on this extra energy 
would not adequately discriminate against signal-like background which is due to partially reconstructed events. 

The overall efficiencies of the basic selection criteria described above are 38% and 30% for the tt" and r? analyses, 
respectively. These estimates have been obtained using MC signal events generated in the range between 1.5 
and 9 GeV^ 
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2. Background Suppression 



The background conditions in two-photon events of low particle multiplicities with tags detected at large and 
(relatively) small polar angles are different. To provide an adequate background suppression for both regions of polar 
angle, wc separate signal event candidates into two samples that have undergone different experimental cuts. In 
this subsection we describe this event separation, the sources of random background and the event selection criteria 
applied to each sample to suppress random background. 

When the scattering angle of the tag is larger than 24.5° (as determined from the calorimeter) we select events 
that have been triggered by the track-based LO trigger. In addition, we require that these events have exactly one 
reconstructed charged track consistent with the tag's shower. There is no efficiency loss associated with the tracking 
requirement which discriminates against background arising mainly from radiative Bhabha events accompanied by 
photon conversion or bremsstrahlung. We include these events in the track-tagged sample. When the scattering angle 
is less than 24.5° we accept both track- and energy-based LO triggers and do not require the presence of the tag's 
track, because the efficiencies of the track-based LO trigger and track reconstruction vanish for tags detected in this 
region of polar angles. Wc include these events that have been triggered either by the track-based or energy-based 
LO trigger in the energy-tagged sample. Notice that while the events from the track-tagged sample must be track- 
triggered, the events from the energy-tagged sample could be either track- or energy-triggered. Tracking information 
for events from the track-tagged sample is utilized in background estimates. The track reconstruction efficiency for 
energy-triggered events is zero. 

Before imposing further selection criteria we obtain improved estimates of the tag energy and direction by using 
transverse-momentum balance and the tag coordinates in the calorimeter. The transverse momenta of the tag and 
of the photon pair should be nearly identical for signal events because the untagged electron usually carries very 
little transverse momentum (below 5 MeV/c) according to the prediction of the MC simulation. Since the transverse 
momentum of the photon pair is measured with much better precision than that of the tag, we equate the magnitude 
of the transverse momentum of the tag with that of the signal photon pair. To calculate the direction of the tag we 
require that its trajectory in the magnetic field goes through the center of the tag's shower. To estimate the center of 
the shower we use the measurement from the calorimeter when the shower is found at polar angles larger than 16.5° . 
At smaller polar angles, however, we use the geometrical center of the crystal with the largest detected energy. This 
is necessary in order to reduce the discrepancies between the data and MC simulation. Using the estimates of the 
tag energy and direction obtained from the transverse-momentum balance we estimate the missing energy and the 
magnitude of missing momentum. To suppress the background from partially reconstructed events we select events 
where the discrepancy between the missing energy and missing momentum is less than 2.3 GeV. This cut is 98% 
efficient for our signal. Wc note that at a later stage of the analysis procedure we will obtain more precise estimates 
of the tag energy and direction. 

Not only should the magnitudes of the transverse momenta of the tag and the photon pair be nearly equal, 
their directions are expected to be practically opposite in the plane perpendicular to the beam collision axis. We 
use the acoplanarity angle, which is the deviation from this expectation, to suppress the background arising from 
radiative Bhabha events with bremsstrahlung photons produced in the materials of the detector. An event of this 
origin enters the energy-tagged sample when the track-based LO trigger is inefficient and a track associated with 
an electron which radiated in the barrel part of the detector cannot be reconstructed. While for signal events the 
acoplanarity distribution peaks near zero, for background events it peaks around 12° for the CLEO geometry and 
CESR kinematics. Acoplanarity discriminates between signal and background events because the measured angular 
position of the shower created by the electron that has undergone bremsstrahlung is shifted with respect to its 
direction at the primary interaction point. This shift is due to the bending of the electron track in the magnetic 
field. To suppress this random QED background in the energy-tagged sample, we select events with acoplanarity less 
than 5°. The background rejection power of this cut exceeds 10, while efficiency loss varies between 20% and 10% for 
between 1.5 and 2.5 GeV^. For larger than 2.5 GeV'^ the efficiency of the acoplanarity cut for the energy-tagged 
events is 90%. In contrast to the energy-tagged sample, the track-tagged sample contains very few bremsstrahlung- 
accompanied radiative Bhabha events because each of these background events has an additional charged track and 
does not pass basic selection criteria (the track reconstruction efficiency for high-energy electrons detected in the 
barrel part of the detector is practically 100% for events recorded by the track-based LO trigger). We select the 
track-tagged events with acoplanarity less than 15°. The efficiency of this loose cut on acoplanarity is 99%. 

We use the decay angle 9d to further suppress background arising from radiative Bhabha events accompanied by 
low-energy split-off clusters. The decay angle is determined from the directions of the 7r° (or 77) candidate in the 
lab frame and one of the daughter photons in the center-of-mass frame of tt" (or rf). Simulation of the detector 
acceptance predicts that the distribution of |cos6'd| is flat between 0.0 to 0.95 and decreases rapidly beyond 0.95 
due to the acceptance loss for soft photons. In contrast to the signal, radiative Bhabha events with split-off clusters 
congregate at | cos^d| = 1.0 because these clusters typically are of low-energy. We reject these asymmetric decays by 
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requiring jcos^d] < 0.90. 

The acoplanarity and decay angle cuts do not eliminate random background completely, because radiative Bhabha 
events accompanied by 7 conversions in detector materials look similar to signal events when triggered by the energy- 
based LO trigger. However, we have found that the shape of this background is monotonic within the signal and 
sideband regions of the 77 invariant mass distribution in both analyses. 

3. Event Quality Requirements 

The angular spectrum of the scattered electrons peaks sharply at small polar angles due to the kinematics of 
processes studied in our analyses. Thus, to measure the cross sections for two-photon production in a tagged mode 
we must understand this critical region of our experimental apparatus very well. While we can, in principle, detect 
tags at polar angles as small as 13° , the fraction of the tag energy collected in the calorimeter at these small polar 
angles is usually less than 20% and might be insufficient to trigger an event. In addition, even if the trigger is satisfied, 
an event might be rejected by the LVL3 filter, which is biased against events with large net transverse momenta. To 
select events identified in the detector regions where the trigger and LVL3 efficiencies are well understood, we need 
better estimates of the tag energy and scattering angle. 

To make precise estimates of the tag energy and scattering angle we use energy-momentum conservation assuming 
that the only particle missing detection is the untagged electron with zero transverse momentum. In practice, this 
method allows us to estimate the parameters of the tag when we measure only the four-momentum of the hadronic 
system and assume that we know the charge of the untagged electron (from crude measurement of the direction of 
missing momentum). From conservation laws we estimate the tag energy E with an r.m.s. resolution of 0.003 E 
and the scattering angle with an r.m.s. resolution of better than 0.6°. In addition, to estimate the scattering angle 
for track-tagged events we use the polar angle of the reconstructed charged track associated with the tag. By using 
the polar angle of the track we achieve an additional small improvement in the resolution of the scattering angle for 
these events. The tag energy for track-tagged events, however, is estimated from energy-momentum conservation; 
i.e. no tracking information is used to estimate the tag energy. In further discussions the values of the tag's 
parameters estimated from energy-momentum conservation and the polar angle of the charged track are referred 
to as constrained values of the tag energy and scattering angle. In Fig. |l| we show the resolution functions of the 
tag energy and scattering angle determined from the differences between analyzed (i.e. measured or constrained) 
and generated quantities (normalized to the generated value for the energy resolution function). These resolution 
functions have been obtained using simulated vr'' events which have values of between 1.5 and 9 GeV^ and satisfy 
all selection criteria discussed above. In our analyses we estimate for each event using constrained values of the 
tag energy and scattering angle. This results in an r.m.s. resolution that varies between 0.1 and 0.3 GeV^ for the 
region between 1.5 and 9 GeV^. 

To isolate the detector region for which the efficiency is small and poorly understood, events with constrained values 
of the tag scattering angle less than 15° are rejected from further analysis. In addition, to reduce the systematic 
uncertainty in the efficiency of the LVL3 filter we select events in which the detected fraction of the tag energy is 
at least 50%. This fraction is estimated from the calorimeter measurement and the constrained value of the tag 
energy. The efficiency of this fractional tag-energy cut is 90% for tags which scatter at 15° and is practically 100% 
for tags which scatter at angles larger than 19°. We have measured the dependence of this efficiency on the polar 
angle using radiative Bhabha events triggered inclusively by the barrel TF-based LO trigger. We show the efficiency 
of the fractional tag-energy cut in Fig. H. 

4. Event Selection Results 

In Figs. ^ and ^ we show the 77 invariant mass distributions for data events that pass all selection criteria for the 
tt" and rj candidates and have values of between 1.5 and 9 GeV^. The points with error bars in these figures 
represent event yields in data. The solid line in each figure shows the result of the binned likelihood fit to data with 
the signal line shape obtained from the MC simulation and an approximation of the remaining random background. 
In the 7r° — > 77 analysis, the background arising from radiative Bhabha events accompanied by photon conversions 
is approximated by an exponential. In the r; — > 77 analysis random background is approximated by the sum of 
an exponential and a constant because the 77-mass distribution shown in Fig. ^ contains two major background 
components. While the first component has the same source as in the tt" analysis, the second component is due to 
radiative Bhabha events with bremsstrahlung radiation in the interface between the drift chambers. 
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C. Background Estimates 



The data may contain tt and rj events that are due to beam-gas interactions or partially reconstructed events 
of higher particle multiplicities. To estimate the beam-gas contribution we use the distributions of the event vertex 
position, visible energy and squared missing mass. Given the profile of the residual gas density near the beam-collision 
point, the vertex position of beam-gas events is much more diffuse than that of the signal. In addition, while beam-gas 
events should have visible energy (i.e., total energy detected in the calorimeter) less than the beam energy, most 
events from the track-tagged sample have visible energy larger than the beam energy. However, at small scattering 
angles the tag needs to go through a larger amount of the detector materials than at large scattering angles and 
can lose a significant part of its energy before reaching the calorimeter. As a result, a large fraction of events from 
the energy-tagged sample (about 20%) falls into the visible energy region below the beam energy. For these energy- 
tagged events we have studied the distribution of the squared missing mass estimated assuming the electroproduction 
hypothesis e^p e^p7r'^( or r]). Using the discriminating power of the distributions described above we conclude 
that the beam-gas background is very small and warrants no subtraction. 

To estimate the background contribution to the track-tagged sample due to e'^e^ annihilation we have studied 
the correlation between the charge and the direction of the tag's track. Signal processes should produce virtu- 
ally all positrons in the +z hemisphere and electrons in —z hemisphere, where +z is the direction of the positron 
beam. However, e"'"e~ annihilation should produce practically the same number of electrons (and positrons) in both 
z-hemispheres. We do not observe a single data event in which this charge-direction correlation indicates e"^e~ an- 
nihilation processes. We conclude that the background from e'^e~ annihilation is fewer than 1 event in both track- 
and energy-tagged samples because the angular distribution of the electrons from this background source is expected 
to be relatively uniform (compared to the rapidly changing signal). 

Finally, there may be some background from other single-tagged two-photon processes. The process 
e^e~ — > e+e~/2(1270) followed by the decay /2(1270) — > tt^tt" is the most likely source of the feed-down for the 
tt'' 77 analysis. To estimate the feed-down from this process, we remove the cut on the energy of the fourth, 
least energetic cluster and repeat the analysis. We estimate that out of 1300 7r° event candidates in data, 80 ± 40 
events are due to the feed-down, where the error refiects the uncertainty of our method. This uncertainty arises 
from the fact that the tt" misidentification probability for the feed-down from the decay /2(1270) depends on 

the relative strengths of the couplings between the tensor meson and two space-like photons of various total helicity 
(determined in the center-of-mass frame of /2(1270)). The central value of the background estimate quoted above 
has been derived assuming that the /2(1270) production proceeds exclusively via the helicity ±2 channel. The error 
refiects the uncertainty in the background estimate which becomes larger (smaller) when we assume that /2(1270) 
is produced only in the helicity (±1) state. We assign this large error to the background estimate because the 
contributions of different helicity amplitudes to the single-tagged cross section for this background process have not 
been measured yet. We observe the /2(1270) feed-down at below 4 GeV^ and subtract its contribution to each 
interval using the shapes of the energy spectra of an additional cluster measured from data and signal MC simula- 
tion. We do not observe a feed-down in the rj —> 'y^f analysis. We have also studied the feed-down from single-tagged 
two-photon processes of higher final-state particle multiplicities such as the production of rj and rj' and estimate the 
overall contribution from these background processes to be insignificant in both analyses. 



D. Systematics 

Contributions to the systematic errors arise from four sources. The primary uncertainty is due to systematic biases 
in the determination of the event selection efficiency. These biases are detailed below. The second contribution is 
a 1% systematic error on integrated luminosity |^ . This error is based on estimates of the theoretical uncertainties 
in the QED radiative corrections in the MC event generators for the processes e^e~ e^e~ and e^e~ — > 77 which 
are employed in the determination of integrated luminosity. The third contribution is a 1% systematic error due 
to the background estimation procedure. The fourth source of systematic error is due to small uncertainties in the 
branching fractions for studied decay chains. This error is negligible in the tt" — > 77 analysis and is less than 1% in 
the — > 77 analysis. 

The largest systematic error is due to the fractional tag-energy cut. We have measured the efficiency of this cut 
using radiative Bhabha events in data. The relative statistical error in this efficiency is less than 3% for polar angles 
larger than 15° so we conservatively include a 3% error to the systematics of energy-tagged events. Note that the 
fractional tag-energy cut is fully efficient for track-tagged events, so no contribution is made to their systematics. 

The efficiency of the LVL3 filter has been measured using tt" signal data events that would have normally been 
discarded by this filter. The statistical error in the measured efficiency is 2% and this gives an estimate of the 
systematic error. 
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The next error comes from the uncertainty in the photon reconstruction efficiency. We have determined this 
uncertainty to be 2%, or 1% per photon from a global fit of the measured ratios of the rj and 77' branching fractions 
to their average values js^. 

We have measured the efficiency of the VD LO trigger over the entire data sample using the TF-triggered endcap 
Bhabha events and have found that this efficiency varies by up to 2.5% of its central value between data subsets. In 
our analysis we use the average value for the VD LO trigger efficiency of 80% and include its r.m.s. variation of 2% 
to the systematic error for track-tagged events. 

To estimate the systematic uncertainty in the efficiency of the extra energy cut we have utilized the shape of the 
extra energy distribution measured from signal data in the rj' — *■ 67 analysis. We estimate this uncertainty to be 2%. 

The efficiency of the acoplanarity cut for energy-tagged events is between 80% and 90%, depending on Q^. To 
estimate the uncertainty in this cut, we have measured its efficiency assuming that the detector simulation system- 
atically underestimates or overestimates azimuthal angular positions of all showers and the tag in the calorimeter by 
one standard deviation of the angular resolution function. We find that under these conditions the efficiency varies 
by less than 1% of itself in any interval. We include this value of 1% to the systematics of energy-tagged events. 

We have also studied other sources of uncertainties such as the efficiencies of missing energy-momentum and decay 
angle cuts and conclude that their total contribution to the systematics is insignificant. 

We include the systematic uncertainties in the amount of feed-down background and in the shape of the 77-mass 
spectrum for random background to the statistical error on the number of signal events in each interval. These 
errors are between 1% and 5% being larger at smaller Q^. 

While the acoplanarity and fractional tag-energy cuts afi^ect only energy-tagged events, the track-based LO trigger 
is specific for the track-tagged events. Thus, the systematic uncertainties associated with the two event samples are 
different. To estimate the systematics for each interval we have used the distributions for MC events which 
belong to the energy- and track-tagged samples. We show these distributions in Fig. ^ 

Our analyses should not be significantly afi'ected by the QED radiative corrections. To order a^, in addition to 
the vacuum polarization and one virtual photon exchange, these corrections describe the processes e'^e~ — > e^e~TZ'y, 
where TZ is one of the studied pseudoscalar mesons |62|. When a radiative photon carries away part of the initial 
center-of-mass energy and remains undetected, we use the nominal value of the beam energy before the radiation 
and overestimate according to Eqn. ^ However, when we estimate the tag energy and scattering angle from 
energy- momentum conservation, we underestimate . Both distortions described above are small effects because 
the energy spectrum of radiative photons is very soft. We neglect the effect of the QED radiative corrections on the 
smearing of the spectrum because these two small effects largely cancel each other. The net smearing is such that 
in our analysis procedure the measured cross sections are insignificantly underestimated. 

There is another aspect of the QED radiative corrections that might need to be taken into account. Namely, 
when we unfold the differential cross sections and obtain the transition form factors, we rely on the prediction of 
a numerical integration that does not contain these corrections and underestimates the cross sections. We expect 
the QED radiative corrections to the cross sections for single-tagged events to be smaller than 2.4% |62|6^ and this 
gives a 1.2% estimate of the systematic uncertainty introduced in the values of T-y'-fU from the unfolding procedure]^. 
Finally, we should emphasize that in order to account for the QED radiative corrections in a consistent manner we 
should have had these corrections implemented in the MC event generator that we use to measure the detection 
efficiency. We did not use such an event generator in our analysis. 

The efficiencies of the event selection criteria employed in our analysis are not fiat over the studied region. 
Most systematic errors for these efficiencies are quoted for a region of low {i.e. less than 3 GeV^) where the 
efficiencies are smaller and the systematic uncertainties are larger than at high (above 3 GeV^). These estimates 
are conservative in the high region where a small fraction of signal events has been detected. 

In the analyses of 77 final states the systematic errors contribute a 5% uncertainty to the measured cross sections. 
As we described above, this uncertainty includes a contribution of ~ 3% that comes from different sources for energy- 
and track-tagged event samples. 



^ If we include the corrections that are due to the vacuum polarization of the probe [i.e. highly virtual) photon in 
the definition of the measured form factors, the remaining QED radiative corrections to these form factors would be 
smaller than 0.5%. The vacuum polarization and all other corrections are of opposite signs and partially cancel each 
other. 
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IV. ANALYSES OF SINGLE- TAGGED 67 AND IO7 FINAL STATES 



The following subsections mainly describe the differences among the analyses of 67 and IO7 final states and the 
previously described analyses of 27 final states, since they share many common features. 



In addition to the trigger utilized for 77 final states, 67 and IO7 single-tagged events have been collected with a 
modified energy-based LI trigger, which is fulfilled when the high energy shower associated with the tag candidate 
is found in the endcap calorimeter and two well-separated clusters, each of detected energy above 100 MeV, are 
identified in the barrel calorimeter. This additional trigger option is especially important for 37r''^67 and 57r°^107 
final states because few of the photons resulting from the vr" decays have sufficient energy to satisfy the high-energy 
trigger threshold of about 500 MeV. 

Each event candidate should contain a tag (in the endcap part of the calorimeter), six, seven, ten or eleven photon 
candidates, and no charged tracks except the tag's track, if reconstructed. The efficiencies of these basic selection 
requirements are 30%, 31% and 12% for the r; —> 67, 77' —> 67, and 77' — > IO7 analyses, respectively, with the 
reconstruction efficiency of about 80% per photon being the dominant source. 

To reduce the systematic uncertainty in the trigger efficiency we select events with the most energetic photon 
candidate detected in the barrel calorimeter at polar angles above 37°. The trigger efficiency for these events is 
larger than 90%. We apply the same missing energy-momentum cut of 2.3 GeV as in the 77 analyses. We require 
acoplanarity less than 30° and do not apply a decay angle cut because there is no need to suppress the small 
background due to radiative Bhabha events. 

Only events that contain at least one combination of the required number of it" 77 and 7; ^ 77 candidates are 
accepted for further analysis. To give an example, we consider the decay chain ri 3tv° 67. Among six or seven 
photon candidates, there must be at least one set of three vr" candidates, where each vr" candidate is identified within 
[-9.0, 3.5](7 of the nominal tt" mass. The mass resolution a has been measured as a function of energy and polar angle 
from data, with a typical value between 6 and 8 MeV/c^. If there is more than one way to form three tt" candidates, 
we use the best combination, i.e. the one which has the smallest x^, where 



We follow the same procedure for 67 and IO7 final states in which we search for the best tt^tt^j; and Stt" combinations, 
respectively. To obtain a better estimate of the parent particle four-momentum we perform a kinematic fit for each 
77-decay candidate from the best combination. For events in which we find an additional energy cluster that has not 
been used to form any of the 7r° or rj candidates, we require that the energy of this cluster be less than 200 MeV. In 
contrast to the 77 analysis, this energy cluster is not necessarily the least energetic one. 

Events that are accepted for further analysis must have constrained values of the tag scattering angle larger 
than 15°. In addition, the detected fraction of the tag energy must be at least 50%. To estimate the constrained 
values of the tag energy and scattering angle we employ energy and momentum conservation laws in which we use 
the four-momenta of the reconstructed tt" — > 77 and 77 — > 77 candidates obtained from the kinematic fits. 

In Figs. ^|-^ we show the invariant mass distributions for data events that pass all selection criteria for the r] and 77' 
candidates and have values of between 1.5 and 9 GeV'^. The points with error bars in these figures represent event 
yields in data. The solid line in each figure shows the result of the binned likelihood fit to data with the signal line 
shape obtained from the MC simulation and a first-order polynomial chosen to approximate the remaining random 
background. 



To estimate the feed-down background, we have studied the distribution of extra energy when the cut on this 
quantity has been removed. We conclude that out of 187 event candidates for the decay 77 3^" in data, 7 events are 
due to feed-down from the decay chain 77' — » 7r°7r°?7 — » 57r°. To subtract this feed-down background, we use the extra 
energy spectra measured from data and signal MC simulation. We do not observe a feed-down in the 77' analyses. 
We estimate the beam-gas and e"*'e~ annihilation backgrounds to be less than 1% of the signal in each analysis. 

In the analyses of 67 and IO7 final states we include a 1% error to the systematics due to the uncertainty in the 
efficiency of the barrel energy-based LI trigger. To estimate this uncertainty we have studied the efficiency of a 
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low-energy trigger threshold for signal data and MC events which have been inclusively triggered with a high-energy 
trigger threshold. All other systematic uncertainties have been discussed in Section III.D. 

In the analyses of 67 and IO7 final states the overall systematic uncertainties in the measured cross sections are 7% 
and 11%, respectively. 



V. ANALYSES OF SINGLE-TAGGED FINAL STATES WITH CHARGED PIONS 

In this section we describe the analyses of final states that contain the tag, two or four charged pions, and at least 
one photon. 



A. Trigger 

As we described in preceding sections, charged tracks can be reconstructed only in events which have been recorded 
with the track-based LO trigger. This trigger is satisfied by two well-separated TF hits, or one TF hit and a VD 
track. The LO triggers are not correlated with the LI triggers; when any of the LO triggers is satisfied, all LI triggers 
are examined js^ . 

In addition to the energy-based LI trigger described previously, there are several track-based LI triggers which are 
efficient for events with charged particles: 

• The LI "electron" trigger is satisfied by a high-threshold bit in the barrel calorimeter and a charged track 
penetrating more than halfway through the volume of the main drift chamber. 

• The LI "two-track" trigger is efficient for events with two or more low transverse-momentum charged particles; 
it requires at least two hits in either region of the TF system, two well-separated low-threshold clusters in the 
barrel calorimeter, and two charged tracks, each of transverse momentum above 90 MeV/c. The L2 trigger is 
fulfilled when at least one charged track of transverse momentum larger than 340 MeV/c is identified. 

• The "hadronic" triggers are designed for multi-particle final states from e^e~ annihilation, but have significant 
efficiency for this analysis as well. These have a variety of possible criteria involving the drift chambers, TF, 
and low-threshold bits of the calorimeter. In general, at least three tracks are required. 

Associated with these track-based LI triggers, earlier data sets had a L2 requirement of a VD hit pattern consistent 
with a charged track of transverse momentum larger than 125 MeV/c. The LVL3 filter does not reject events that 
are collected by the track-based LI triggers. 



B. Analysis Procedure 

Each event candidate must contain the tag, an exact number of charged tracks (excluding the tag's track, if 
reconstructed), and at least as many photon candidates as are needed for full reconstruction of a studied decay chain. 
All tracks except for the tag's track are assumed to be due to charged pions. The net charge of the reconstructed pions 
must be zero. Photon candidates include all barrel (endcap) calorimeter clusters of energies larger than 30 (50) MeV 
except for those that are closest to the intersection points of charged tracks with the calorimeter. The efficiencies of 
the basic requirements described above are defined by the charged pion and photon reconstruction efficiencies, each 
about 80% per particle. 

To select events that trigger with high efficiency and small systematic uncertainty, we impose several event-quality 
criteria. Namely, we require that at least one charged track of transverse momentum larger than 250 MeV/c be 
detected. In addition, the charged track of largest transverse momentum and either the tag or the most energetic 
photon candidate must be detected in the barrel calorimeter at polar angles above 37°. Finally, we reject events 
which contain charged tracks of momenta less than 80 MeV/c because for these tracks the systematic uncertainty in 
the track reconstruction efficiency is large. 

Given that p°7 events are primarily recorded with the energy-based LI trigger, tighter event selection criteria 
are imposed in this analysis. We select events which have at least one charged track of transverse momentum 
above 450 MeV/c. The most energetic photon candidate must have energy, E-y, larger than 130 MeV. We assume 
this photon candidate to be due to the signal process ri' — > p'-^'y. To suppress random background we select events 
with the reconstructed -k'^tv~ mass between 550 and 800 MeV/c^. This is referred to as the p^-mass cut. 

In the analyses of the final states that contain the decays 7r° — > 77 and rj — > 77, events must contain at least one 
combination of the exact number of the candidates for these decays as required for full reconstruction of the studied 
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decay chain. The energy clusters that enter the best combination are assumed to be signal photons. The total energy 
collected in the calorimeter clusters other than the signal photon candidates and the energy clusters matched to the 
projections of the charged tracks must be less than 500 MeV. These extra energy clusters are mostly due to the 
interactions of the charged pious with the materials of the detector. No requirement is made on the number of such 
clusters. 

We use the momenta of the charged tracks and signal photons (after kinematic fits, where applicable) and employ 
energy-momentum conservation to estimate the tag energy and scattering angle. We select events in which the 
detected fraction of the tag energy is at least 50% and a scattering angle is larger than 15° where both parameters 
are estimated using energy-momentum conservation. 

In the T] — > ■iv'^n~-K° analysis we need to suppress a large feed-down from the decay chain rj' ■k'''it'^t] — > 7r^7r~37r°. 
To suppress this feed-down, we require the difference between the measured and constrained values of the tag scattering 
angle be less than 2°. The feed-down suppression power of the combination of this and the extra energy cuts is a 
factor of 23, while the efficiency loss is less than 3%. 

We employ the particle identification capabilities of our apparatus to reduce the large random background observed 
in the rj' — > p°7 analysis. This random background is primarily due to the process e'^e~ — > e^e~e^e~ accompanied 
by bremsstrahlung radiation, split-off showers or beam-related energy clusters. To suppress random background 
we utilize the fact that specific ionization energy losses, dE/dx, are larger for electrons than for charged pious. 
This information is used in the requirement on -Pj^2, the upper tail probability of the distribution of the dE/dx 
measurements for charged pion candidates In the ideal case (i.e. if the dE/dx distribution were Gaussian) the 
correct choice of the particle-identification hypothesis would produce a uniform P^2 distribution, while events with 
an incorrect particle-identification hypothesis tend to congregate near zero. We calculate P^2 for the tracks assuming 
them to be due to charged pious. To suppress unwanted background events, P^2 is required to be larger than 0.005. 
The efficiency of the P^2 cut is not 99.5% but 98% because a small fraction of the signal events (in both data and 
simulation) does not have dE/dx information and the dE/dx distribution has non-Gaussian tails. The same cut on 
P^2 is applied in all analyses with charged pious. 

In Figs, [gj-p^ we show the invariant mass distributions for data events that pass all selection criteria for the r) 
{rj') candidates and have values of between 1.5 and 20 (30) GeV^. The points with error bars in these figures 
represent event yields in data. The solid line in each figure shows the result of the binned likelihood fit to data 
with the signal line shape obtained from the MC simulation and a linear approximation of the remaining random 
background. The remaining random background observed in the analysis of the 7r^7r~7 final state is due to the 
process e^e~ e^e~ ji^ ii~ accompanied by noise and split-off energy clusters. 



C. Background Estimates and Systematics 

To estimate the feed-down background, we have analyzed the distributions of extra energy and the difference 
between the measured and constrained values of the tag scattering angle when the cuts on these quantities have been 
removed. We conclude that fewer than 2 events in the r) 7r"'"7r~7r'' analysis are due to feed-down from the decay 
chain rj' ■K^n'^ri 7r^7r~37r°. We have not identified any feed-down background in the rj' analyses. We estimate 
that the background contribution from beam-gas interactions and e'*'e~ annihilation processes is less than 1% of the 
signal in all analyses. This gives an estimate of the relevant systematic uncertainty. 

To estimate the systematic uncertainty in the efficiency of the LO trigger, we select signal events that are triggered 
by the TF-based LO trigger and measure the VD efficiency per event. Using a similar method we measure the efficiency 
of the TF-based LO trigger for events that are triggered by the track-based LO trigger. We estimate the uncertainty 
in the LO trigger efficiency to be 1%, which is the typical deviation between either of these efficiencies measured from 
data and simulation. Note that the efficiency of the LO trigger is higher than 98% for events which satisfy the basic 
selection requirements. 

We have measured the efficiency of the P^2 cut in a nearly background-free environment using fitted mass distribu- 
tions for signal events in data and MC simulation for the decay chain rj' — ^ -K^n'Tj — > 7r^7r~27. We have found this 
efficiency to be 98% (as discussed in the previous subsection) and use the 2% statistical error of this measurement as 
an estimate of the systematic uncertainty. 

To estimate the uncertainty in the efficiency of the E-y cut in the rj' p^j analysis, we have measured this efficiency 
assuming that the energies of the reconstructed photons in the simulation are systematically shifted by 2% of their 
nominal values. We have observed a relative change of 1% in the efficiency of the E-y cut and this gives an estimate 
of its systematic uncertainty. 

We estimate the uncertainty in the track reconstruction efficiency to be 2% per charged pion. It is determined 
from a global fit of the measured ratios of the 77 and rj' branching fractions to their average values kJ. 
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The uncertainty in the efficiency of the p'^-mass cut is negligible because, except for the p^-line shape, the matrix 
element for the decay chain r( ^"7 7r^7r~7 is determined by QED and kinematics. To confirm this statement 
we remove the p''-mass and cuts and compare the distributions of E-^ and | cos Q* \ measured from signal data 
and MC simulation, where E-^ is the signal photon energy in the lab frame^. These distributions are shown in 
Figs. |l| and ^. We observe good agreement between the data and MC spectra of E^ and | cos Q* \ and conclude that 
the approximations given by Eqns. ^ and ^ describe the data well. We note that both figures show the observed 
spectra, i. e. no detection efficiency corrections have been applied to these distributions. The good agreement between 
the shape of the | cos 6* \ distribution obtained from the simulation and sin^ 6* curve is due to the detection efficiency 
being practically flat over the full range of | cos 9* \ . 

All other systematic uncertainties have been discussed in Sections III.D and IV. B. In the analyses of final states with 
charged pions the overall systematic uncertainty in the measured cross sections is between 7% and 10%, depending 
on the final state. 



VI. UNFOLDING PROCEDURE FOR THE TRANSITION FORM FACTORS 

To measure the products of the differential cross sections and branching fractions for each decay chain we use 
the following analysis procedure. Data events that pass all selection criteria are used to form the distribution 
where the value of for each event is estimated from energy-momentum conservation (and the polar angle of the 
tag's track when the track is reconstructed) . Next we divide the event yields into intervals. For each interval 
we obtain the number of signal events in data from the fit to the invariant mass distribution. Then we estimate 
and subtract the feed-down background using the methods described in preceding sections. Finally we correct the 
background-subtracted number of signal events for the detection efficiency. The signal line shapes used in the fits 
and the detection efficiencies are determined from the detector simulation for each interval. 

To extract the transition form factors we compare the measured and the predicted values of the cross sections. 
Namely, for each interval, we measure the form factors J-tli*""^ (0^) from: 

\T',i^\&t - ^^\^''UQ')\^ (12) 

where J-^l^'fi{Q^) is the approximation for the Q'^-dependent part of the form factor in MC simulation, and a{data) 
and a{MC) are the cross sections for this interval measured in data and predicted using numerical integration, 
respectively. The transition form factors are measured at where the differential cross sections achieve their mean 
values according to the results of numerical integration. The numerical results have been obtained at an average 
center-of-mass energy of 10.56 GeV with the approximation for the form factor given by Eqn. ^ 

The distributions measured from data and obtained numerically are shown in Figs. ^ and |l7| for the tt" — > 77 
and rj — > 77 analyses, respectively. Only statistical errors are shown in these figures. To plot the results of numerical 
integration we use r(7r° — > 77) = 7.74 eV and T{ri — > 77) = 463 eV |59) . 

We show our experimental results in Tables I-X. These tables show the intervals, event yields obtained from the 
fits, numbers of signal events after subtraction of the feed-down background, detection efficiencies, the values, the 
products of the differential cross sections and relevant branching fractions, and the transition form factors, represented 
in the form Q^\J-'y»'j-R.{Q^)\. In Tables I-X the first error is statistical and the second error (where given) is systematic. 

VII. COMPARISON OF THE RESULTS WITH THEORETICAL PREDICTIONS 

In this section we compare the results for tt*^ with theoretical predictions. For the transition form factors of and 
rj' we compare the results with the PQCD asymptotic prediction only because little is known in theory about the 
wave functions of these mesons. No predictions for the form factors of rj and r/' are available at this time except 
for the prediction of KroU et al. ||l^ where these authors assumed that the shapes of the wave functions of all three 
pseudoscalar mesons are similar. 



These figures of merit for the analysis of the decay chain 77' — » p''7 tt+tt 7 were proposed in 
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A. Results for n 

In Figures we compare our results for Q^I^*^^o(Q^)| with the theoretical predictions. Also shown in these 

figures are the results of the CELLO experiment jHoj and the asymptotic prediction of PQCD given by Eqn. ^. For 
both experimental results the error bars represent the statistical errors only. To plot the results of the theoretical 
predictions we use their published analytical forms. To estimate the value of we use Eqns. I and § and the 
tabulated two-photon partial width of vr" [^. This estimate of (92.3 MeV) agrees with its experimental value 
(92.4 MeV) which has been measured previously from charged pion decayij^ pg[ ]. 

In Fig. ^ the results are compared with the predictions made by Jakob et al. . These authors calculated the 
7*7 —> 7?. transition form factor by employing a PQCD-based technique and QCD radiative corrections [^. They 
used two estimates for the 7r° wave function: the asymptotic wave function and the CZ wave function. This theoretical 
prediction gives a much better agreement with our results when the asymptotic wave function is used. In terms of 
the PQCD-based approach this indicates that the wave function has already evolved to the asymptotic form at as 
small as 1 GeV'^. Notice that T^-^^o calculated with the CZ wave function changes when the QCD evolution of this 
wave function over the studied range is taken into account according to The transition form factor does not 
change when the asymptotic wave function is used because this wave function exhibits no QCD evolution to leading 
order in Ua- However, in next-to-leading order in Os any wave function, including the asymptotic, is subject to the 
QCD evolution [^. If this evolution is taken into account the prediction with the asymptotic wave function which 
has been derived to leading order in would also change slightly prf . 

Cao et al. also made a prediction based on PQCD . These authors disagreed with the approximations made to 
simplify the form of HSA in Their prediction includes transverse momentum corrections and is compared with 
our results in Fig. ^ for the asymptotic and CZ wave functions. The theoretical prediction of Cao et al. yields a 
smaller value of T^t.^^o for less than 8 GeV^ when the CZ wave function is used. This is a most intriguing result 
because the CZ wave function has been proposed to account for measured excesses in the rates for various processes, 
thus leading to larger values of the form factors and cross sections [^. 

The prediction of Radyushkin et al. based on the QCD sum- rules method is compared with the experimental 
results in Fig. ^ This calculation describes the saturating behavior of our measurement, though it disagrees with 
the data at smaller . It should be noted that at low the prediction is not expected to agree with the data: the 
QCD radiative corrections which would be larger at smaller have not been included in this theoretical analysis. 
The discrepancy between the absolute values of the asymptotic limits of PQCD and of this prediction might be due 
to the uncertainties in the expectation values of the vacuum condensates that are known only with 30% precision . 
However, according to the authors, the agreement can be achieved by means of complicated QCD-evolution analysis 
of the correlator functions used in this theoretical approach |Q. 

Finally, we derive the value of the pole-mass parameter A^o which we use to represent our results in a simple 
phenomenological form. We fit our results for |jF^*^^o((3^)P with a function given by Eqn. ^ and obtain the following 
result: 

A^o = 776 ± 10 ± 12 ± 16 MeV, (13) 

where the first error is statistical, the second error represents systematic uncertainties of our measurements, and the 
third error is due to the uncertainty in the value of r(7r'' 77) The result of the fit is shown in Fig. |2l| While 
we observe that a simple VMD-like approximation describes the data very well, we should note that it disagrees with 
the asymptotic prediction of PQCD. Also shown in Fig. |2^ is the interpolation given by Eqn. |^ 

B. Results for rj 

We show the results of our measurements for Q'^\J--^' -iri{Q'^)\ in Fig. . This figure also shows the asymptotic 
prediction of PQCD given by Eqn. ^ and the interpolation given by Eqn.£ To estimate the value of (97.5 MeV) 
we use Eqns. ^ and ^ and the tabulated two-photon partial width of 77 [£9|. We fit the |.7-"-,*-y,,(Q^)|'^ distributions 
measured using each decay chain with the functional form given by Equ!^ and obtain the values of the pole-mass 
parameter A,, that are shown in Table In this table, for each measurement, the first error is statistical, the 
second error represents systematic uncertainties of our measurement, and the third error refiects the uncertainty in 



For each meson TZ, where 7?. is tt , t; or t;', our definition of the meson decay constant f-n difi^ers by a factor of 
1 / y/2 from the one accepted by the Particle Data Group and given in M . 
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the two-photon partial width of rj. From a simultaneous fit to our three measurements for the production of r; we 
obtain the following value of the pole-mass parameter: 

= 774 ± 11 ± 16 ± 22 MeV. (14) 

The result of this fit is shown in Fig. p2[ 

We use the measured values of the parameters A^o and A^ to compare the soft non-perturbative properties of tt" 
and rj. This is a legitimate comparison because the chiral limit given by Eqn. ^ and the asymptotic prediction given 
by Eqn. ^ are expected to hold for both tt" and rj. From the comparison between the measured values of A^o and A,, 
we conclude that the shapes of the 7*7 meson transition form factors of n" and r; are nearly identical, which 
strongly indicates the similarity between the wave functions of these mesons. 



C. Results for rj' 

We show the results of our measurements for Q^\J--y'-yr]i iQ'^)\ in Fig. ^ This figure also shows what would be 
the PQCD asymptotic prediction given by Eqn. ^ for Q^|i^-y*-™'(Q^)| if the chiral limit given by Eqn. ^ held for rj' . 
To estimate the value of /^/ (74.4 MeV) we use Eqns. W and M and the tabulated two-photon partial width of rj' of 
4.3 keV H. 

We fit the \^j*jr]'{Q^)\'^ distributions measured using each decay chain with the functional form given by Eqn. |^ 

From a simultaneous fit to our six 
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and obtain the values of the pole-mass parameter A^/ that are shown in Table 
results for the production of ri' we obtain the following value of the pole-mass parameter: 

A^/ = 859 ± 9 ± 18 ± 20 MeV. (15) 

The result of this fit is shown in Fig. ^ 

The results of our measurements for the production of rj' demonstrate that if this particle were a qq bound state 
and the QCD chiral limit given by Eqn. M held for this meson, the Q^-dependence of the transition form factor 
of rj' and consequently its wave function would be significantly different from these non-perturbative properties of 
either vr" or 77. 



VIII. CONCLUSIONS 



We have measured the form factors associated with the electromagnetic transitions 7*7 ^ meson in the regions 
of momentum transfer from 1.5 to 9, 20, and 30 GeV^ for the vr", 77, and 77' mesons, respectively. These are the first 
measurements above 2.7 GeV^ for tt" and above 7 GeV^ for 77 and 77'. 

Our measurement for tt" unambiguously distinguishes among various theoretical predictions for the form factors 
of the 7*7 tt" transition. We have demonstrated that the non-perturbative properties of tt" and 77 agree with each 
other which indicates that the wave functions of these two mesons are similar. In the 77' analysis we have shown that 
the non-perturbative properties of 77' differ substantially from those of tt" and 77. Our measurement for 77' provides 
important information for future theoretical investigations of the structure of this particle. 



IX. ACKNOWLEDGEMENTS 



We gratefully acknowledge the effort of the CESR staff in providing us with excellent luminosity and running 
conditions. J.P.A., J.R.P., and I.P.J. S. thank the NYI program of the NSF, M.S. thanks the PFF program of the 
NSF, G.E. thanks the Heisenberg Foundation, K.K.G., M.S., H.N.N. , T.S., and H.Y. thank the OJI program of DOE, 
J.R.P., K.H., M.S. and V.S. thank the A.P. Sloan Foundation, A.W. and R.W. thank the Alexander von Humboldt 
Stiftung, and M.S. thanks Research Corporation for support. This work was supported by the National Science 
Foundation, the U.S. Department of Energy, and the Natural Sciences and Engineering Research Council of Canada. 



[1] V. Savinov, PhD thesis. University of Minnesota, Minneapolis, 1997, Unpublished. 



18 



[2] V.M. Budnev et al, Phys. Rep. C15 (1975) 181. 

[3] G.P. Lepage and S.J. Brodsky, Phys. Rev. D22 (1980) 2157. 

[4] S.J. Brodsky and G.P. Lepage, Phys. Rev. D24 (1981) 1808. 

[5] A.V. Efremov and A.V. Radyushkin, Phys. Lett. B94 (1980) 245. 

[6] V.L. Chernyak and A.R. Zhitnitsky, Phys. Rep. 112 (1984) 173. 

[7] S. Adler, Phys. Rev. 177 (1969) 2426. 

[8] J. Bell and R. Jackiw, Nuovo. Cim. A60 (1969) 47. 

[9] L.H. Ryder, Quantum Field Theory, p. 297, Cambridge University Press, 1986. 
[10] S. Coleman, Aspects of Symmetry, Cambridge University Press, 1985. 
[11] M.A. Shifman, A.I. Vainshtein, V.I. Zakharov, Nucl. Phys. B147 (1979) 385. 
[12] P. KroU, Universitat Wuppertal preprint WU B 94-17. Also, hep-ph 9409262, 1994. 
[13] R. Jakob et al, J. Phys. G22 (1996) 45. 
[14] H.-N. Li and G. Sterman, Nucl. Phys. B381 (1992) 129. 

[15] V. Savinov, CLEO Coll., in the Proceedings of the Conference Photon'95, p. 203, edited by David J. Miller et 

al, World Scientific, 1995. 
[16] F.-G. Cao et al, Phys. Rev. D53 (1996) 6582. 
[17] H.-N. Li, Phys. Rev. D48 (1993) 4243. 
[18] J. Bolz et al, Phys. Lett. B342 (1995) 345. 
[19] M. Dahm et al, Z. Phys. C68 (1995) 595. 
[20] J. Bolz et al, Z. Phys. C66 (1995) 267. 

[21] B. Chibisov and A.R. Zhitnitsky, Phys. Rev. D52 (1995) 5273. 

[22] N. Isgur and C.H. Llewellyn Smith, Phys. Rev. Lett. 52 (1984) 1080. 

[23] J. Botts and G. Sterman, Nucl Phys. B325 (1984) 62. 

[24] N. Isgur and C.H. Llewellyn Smith, Nucl Phys. B317 (1989) 526. 

[25] R. Jakob and P. KroU, Phys. Lett. B315 (1993) 463. 
[26] V. Braun and I. Halpcrin, Phys. Lett. B328 (1994) 457. 

[27] A.V. Radyushkin, in the Proceedings of the Workshop on CEBAF at Higher Energies, p. 273, edited by N. Isgur 

and P. Stoler, CEBAF, 1994. 
[28] A.R. Zhitnitsky, hep-ph 9605226, 1996. 

[29] A.V. Efremov and A.V. Radyushkin, Theor. Math. Phys. 42 (1980) 97. 
[30] A.V. Radyushkin, Nucl Phys. A532 (1991) 141. 

[31] I.V. Musatov and A.V. Radyushkin, Thomas Jefferson National Accelerator Facility (formerly CEBAF) preprint 

JLAB-THY-97-07. Also, hep-ph 9702443, 1997. 
[32] A.V. Radyushkin and R. Ruskov, Phys. Lett. B374 (1996) 173. 
[33] Private communication with A.V. Radyushkin and R. Ruskov. 
[34] P. Ball, Phys. Rev. D48 (1993) 3190. 

[35] A. Khodjamirian and R. Riickl, Preprint WUE-ITP-96-020, also hep-ph 9610367, 1996. 
[36] P. Ball and V.M. Braun, Phys. Rev. D55 (1997) 5561. 
[37] P. Kroll and M. Raulfs, Phys. Lett. B387 (1996) 848. 
[38] V.V. Anisovich et al, Phys. Rev. D55 (1997) 2918. 

[39] A.E. Dorokhov, Universita di Pisa and INFN preprint IFUP - TH 36/95, July 1995. 
[40] S. Ong, Phys. Rev. D52 (1995) 3111. 

[41] R.M. Davidson, to be published in the Proceedings of the Workshop on the Structure of the rj' Meson, edited by 

N. Isgur and P. Stoler, New Mexico State University and CEBAF, Las Cruces, New Mexico, 1996. 
[42] M.R. Frank et al, Phys. Lett. B359 (1995) 17. 
[43] D. Kekez and D. Klabucar, Phys. Lett. B387 (1996) 14. 

[44] A. Anselm et al. Harvard University preprint HUTP-95/A037. Also, hep-ph 9603444, 1996. 
[45] V.M. Bclyacv and M.B. Johnson, hep-ph 9605279, 1996. 

[46] G. Sterman and P. Stoler, Hadronic Form Factors and Perturbative QCD. To be published in Ann. Rev. of Nucl 

and Part. Sc., 1997. 
[47] R.I. Djhelyadin et al, Lepton-G Exp., Phys. Lett. B94 (1980) 548. 
[48] Ch. Berger et al, PLUTO Coll., Phys. Lett. 142 (1984) 125. 
[49] H. Aihara et al, TPC/27 Coll., Phys. Rev. Lett. 64 (1990) 172. 
[50] H.-J. Behrend et al, CELLO Coll., Z. Phys. C49 (1991) 401. 
[51] Y. Kubota et al, CLEO Coll., Nucl Inst, and Meth. A320 (1992) 66. 
[52] C. Bebek et al, CLEO Coll., Nucl Inst, and Meth. A302 (1992) 261. 
[53] D.E. Acosta, PhD thesis. University of California, San Diego, 1993, Unpublished. 
[54] B.H. Ong, PhD thesis. University of California, San Diego, 1993, Unpublished. 

[55] TWOGAM. The Two-Photon Monte Carlo Simulation Program. Written by D.M. Coffman. UnpubUshed. 



19 



[56] J.J. Sakurai, Ann. Phys. 11 (1960) 1. 

[57] A.V. Radyushkin and R. Ruskov, Nucl. Phys. B481 (1996) 625. 
[58] V.N. Gribov and I.Ya. Pomoranchuk, Phys. Rev. Lett. 8 (1962) 343. 

[59] Particle Data Group, L. Montanet et al., Review of Particle Properties, Phys. Rev. D54 (1996) 1. 

[60] "GEANT. Detector Description and Simulation Tool." , by Application Software Group, Computing and Networks 

Division. The CERN Program Library index W5013, CERN, Geneva, Switzerland. 
[61] G. Crawford et al.., CLEG Coll., Nucl. Inst, and Meth. A345 (1994) 429. 
[62] W.L. van Neerven and J.A.M. Vermaseren, Nucl. Phys. B238 (1984) 73. 
[63] M. Landr0 et al, Phys. Rev. D36 (1987) 44. 

[64] M. Abramowitz and LA. Stegun (eds.), Handbook of mathematical functions (National Bureau of Standards, 
Washington; D.C., 1964) The upper tail probability of the distribution (-P^a) is the area of the distribution 
with N degrees of freedom above measured value. In our case the number of degrees of freedom is the number 
of available dE/dx measurements. The CERN Program Library index GlOO, CERN, Geneva, Switzerland. 

[65] D. Miiller, Phys. Rev. D51 (1995) 3855. 



20 



TABLE I. The results of the tt" — > 77 analysis assuming B = £(77° — » 77) = 0.99. The differential cross section is 
for e"'"e~ — > e'''e~7r°. 
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TABLE XL The results of the 77 ^ 77 analysis assuming B = B{ri —^ 77) = 0.39. The differential cross section is 
for e"'"e~ — » e"'"e~7j. 
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TABLE III. The results of the rj 37r° analysis assuming B = B{rj ^ 3n°)x B^{-k° 77) = 0.31. The differential 
cross section is for e^e~ — > e^e~r]. 
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TABLE IV. The results of the r; ^ tv+tt--iv° analysis assuming B = B{r} ^ 7r+7r"7r°)x B{-n° -» 77) = 0.23. The 
differential cross section is for e^e~ — > e^e~ri. 
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21 


235 ± 35 


2.23 


714 ± 105 ± 55 


12.1 ± 0.89 ± 0.47 


2.5 - 3.5 


49 ± 7 


48 ± 7 


23 


210 ± 31 


2.94 


318 ± 47 ± 25 


13.0 ± 0.96 ± 0.50 


3.5 - 5.0 


31 ± 6 


31 ± 6 


26 


117 ± 23 


4.16 


118 ± 23 ± 9 


14.4 ± 1.39 ± 0.55 


5.0 9.0 


32 ± 6 


32 ± 6 


26 


122 ± 23 


6.56 


46 ± 9 ± 4 


20.1 ± 1.88 ± 0.77 


9.0 20.0 


() = :i 


(i = :5 


2.") 


2:! = 10 


12.71 


:5.1 = 1.1 = 0.2 


18.1 = 1.19 = 0.71 



TABLE V. The results of the rj — > tP-k't] 67 analysis assuming B = B{r]' — » tt^tt"?;) x B(r] — » 77) x 
jB^(7r° — > 77) = 0.080. The differential cross section is for e"'"e~ — » e'^e~ri'. 



interval 

(GeV^) 


detected 


signal 


e 

% 


B X N^, 

produced 


Q' 
(GeV^) 


(fb / GeV^) 


(0.01 X GoV) 


1.5 - 2.0 


40 ± 7 


40 ± 7 


8.4 


474 ± 85 


1.73 


4132 ± 740 ± 310 


20.1 ± 1.8 ± 0.8 


2.0 - 2.5 


40 ± 7 


40 ± 7 


16 


259 ± 44 


2.23 


2258 ± 381 ± 169 


22.7 ± 1.9 ± 0.9 


2.5 - 3.5 


29 ± 6 


29 ± 6 


16 


176 ± 38 


2.94 


767 ± 164 ± 58 


21.1 ± 2.3 ± 0.8 


3.5 5.0 


17 ± 4 


17 ± 4 


18 


94 ± 24 


4.16 


274 ± 70 ± 21 


22.7 ± 2.9 ± 0.9 


5.0 - 9.0 


14 ± 4 


14 ± 4 


16 


90 ± 24 


6.56 


98 ± 26 ± 7 


30.0 ± 4.0 ± 1.1 



TABLE VL The results of the 77' n°Tr°ri 5n° IO7 analysis assuming B = B{ri' 7r°7r°r?)x B{ri 3n°)x 
B^{'K° —^ 77) = 0.063. The differential cross section is for e^e~ — > e^e~r]' . 



interval 

(GeV^) 


detected 


signal 


e 

% 


B X N^, 

produced 


(GcV^) 


da/dQ^iQ"") 

(fb / GcV^) 


(0.01 X GeV) 


1.5 - 3.0 


18 ± 5 


18 ± 5 


3.5 


510 ± 153 


2.09 


1875 ± 563 ± 204 


18.6 ± 2.8 ± 1.0 


3.0 - 9.0 


7 ± 3 


7 ± 2 


5.4 


129 ± 49 


4.92 


118 ± 45 ± 13 


20.0 ± 3.8 ± 1.1 



TABLE VIL The results of the r?' -> p°7 7r+7r~7 analysis assuming B = B{r]' -» p°7)x B{p° tt+tt") = 0.30. 
The differential cross section is for e+e" — > e^e~r]' . 



interval 

(GeV^) 


detected 


signal 


e 

% 


B X N^, 

produced 


0" 
(GeV^) 


dGldQ'^{Q'') 
(fb / GeV^) 


(0.01 X GeV) 


1.5 - 2.0 


111 ± 13 


111 ± 13 


8.9 


1257 ± 152 


1.73 


2891 ± 350 ± 197 


16.8 ± 1.02 ± 0.57 


2.0 - 2.5 


131 ± 14 


131 ± 14 


17 


765 ± 84 


2.23 


1759 ± 193 ± 120 


20.0 ± 1.10 ± 0.68 


2.5 - 3.5 


123 ± 14 


123 ± 14 


21 


593 ± 69 


2.94 


681 ± 79 ± 46 


19.9 ± 1.15 ± 0.68 


3.5 - 5.0 


86 ± 11 


86 ± 11 


24 


353 ± 47 


4.16 


270 ± 36 ± 18 


22.6 ± 1.51 ± 0.77 


5.0 - 9.0 


49 ± 10 


49 ± 10 


31 


158 ± 32 


6.56 


45 ± 9 ± 3 


20.4 ± 2.08 ± 0.69 


9.0 - 30.0 


22 ± 8 


22 ± 8 


37 


58 ± 21 


15.30 


3.2 ± 1.1 ± 0.2 


24.8 ± 4.44 ± 0.84 
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TABLE VIII. The results of the r\ tv^-k rj — > -k^tt 27 analysis assuming B = B{r]' — > rj-K'^iv )x 

B{ri —^ 77) = 0.17. The differential cross section is for e~^e~ —^ e^e^r/. 



interval 

(GeV^) 


detected 


signal 


e 

% 


B X N^, 

produced 


(GeV^) 


da/dQ2(Q2) 
(fb / GeV^) 


(0.01 X GeV) 


1.5 - 2.0 


57 ± 8 


57 ± 8 


7.6 


743 ± 104 


1.73 


3007 ± 421 ± 199 


17.1 ± 1.2 ± 0.6 


2.0 - 2.5 


70 ± 9 


70 ± 9 


17 


408 ± 51 


2.23 


1651 ± 208 ± 109 


19.4 ± 1.2 ± 0.6 


2.5 - 3.5 


60 ± 8 


60 ± 8 


21 


282 ± 38 


2.94 


570 ± 77 ± 38 


18.2 ± 1.2 ± 0.6 


3.5 - 5.0 


58 ± 8 


58 ± 8 


27 


216 ± 30 


4.16 


292 ± 40 ± 19 


23.5 ± 1.6 ± 0.8 


5.0 - 9.0 


45 ± 7 


45 ± 7 


34 


133 ± 20 


6.56 


67 ± 10 ± 4 


24.9 ± 1.9 ± 0.8 


9.0 - 30.0 


16 ± 4 


16 ± 4 


36 


44 ± 11 


15.30 


4.3 ± 1.1 ± 0.3 


28.8 ± 3.6 ± 0.9 



TABLE IX. The results of the 7]' — » tt+tt rj — » 27r+27r ■k° — » 27r+27r 27 analysis assuming B = B{r)' — » tt+tt 7?)x 
B{ri — » 7r"'"7r~7r°)x iB(7r° — > 77) = 0.10. The differential cross section is for e^e~ — > e^e~ri' . 



interval 

(GeV^) 


detected 


signal 


e 

% 


B X iV^, 

produced 


(GeV^) 


(fb / GeV^) 


(0.01 X GeV) 


1.5 - 2.5 


33 ± 6 


33 ± 6 


6.2 


528 ± 95 


1.92 


1830 ± 329 ± 176 


15.9 ± 1.4 ± 0.7 


2.5 - 3.5 


22 ± 5 


22 ± 5 


13 


169 ± 40 


2.94 


584 ± 138 ± 56 


18.4 ± 2.2 ± 0.9 


3.5 - 5.0 


18 ± 5 


18 ± 5 


16 


113 ± 30 


4.16 


261 ± 69 ± 25 


22.2 ± 2.9 ± 1.1 


5.0 9.0 


15 ± 4 


15 ± 4 


21 


74 ± 20 


6.56 


64 ± 17 ± 6 


24.4 ± 3.2 ± 1.2 


9.0 - 30.0 


4 ± 2 


4 ± 2 


24 


16 ± 8 


15.30 


2.7 ± 1.4 ± 0.3 


22.9 ± 5.8 ± 1.1 



TABLE X. The results of the r( — > v'^'n 37r — > tt^tt 67 analysis assuming B = {B{ri' —^ n'^n' 
B{r)' — » 7r°7i-°77)x B{r] — > 7r+7r~7r°)) x B^{Tr° — > 77) = 0.14. The differential cross section is for e+e~ — > e+e~r?'. 



interval 

(GeV^) 


Nr,' 

detected 


signal 


e 

% 


B X N„, 

produced 


Q' 

(GeV^) 


rfa/dQ^(Q2) 
(fl) / GcV^) 


Q \J~ -yj^' (Q ) 1 

(0.01 X GoV) 


1.5 - 2.5 


54 ± 8 


54 ± 8 


3.7 


1468 ± 206 


1.92 


2803 ± 393 ± 247 


19.7 ± 1.4 ± 0.9 


2.5 - 3.5 


25 ± 6 


25 ± 6 


7.5 


330 ± 78 


2.94 


630 ± 149 ± 55 


19.1 ± 2.3 ± 0.8 


3.5 - 5.0 


15 ± 4 


15 ± 4 


10 


161 ± 45 


4.16 


205 ± 57 ± 18 


19.7 ± 2.7 ± 0.9 


5.0 - 9.0 


13 ± 4 


13 ± 4 


13 


101 ± 34 


6.56 


48 ± 16 ± 4 


21.0 ± 3.5 ± 0.9 


9.0 - 30.0 


2 ± 1 


2 ± 1 


15 


14 ± 10 


15.30 


1.3 ± 0.9 ± 0.1 


15.7 ± 5.6 ± 0.7 
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TABLE XI. Values of the pole-mass parameters A^o, A^ and A^/ measured using various final states. For each 
measurement, the first error is statistical, the second error represents the systematic uncertainties of our measurement 
and the third error reflects the experimental error in the value of the two-photon partial width of the meson. 



Decay chain 


A,^ (MeV) 


7r" — > 77 


776 ± 10 ± 12 ± 16 


?7 -» 77 


778 ± 19 ± 12 ± 22 


77 — » OTT — > 07 


'7'70 _l_ on _l_ 1 T _l_ 00 
1 16 ± JX) ± 11 ±2,1 




773 ± 18 ± 18 ± 22 


Simultaneous fit to all r; data 


774 ± 11 ± 16 ± 22 


77' — > p"7 — > 7r"'"7r~7 


857 ± 15 ± 19 ± 19 


77' ■K'^TT~ri 7r^7r~27 


864 ± 16 ± 18 ± 19 


77 n n 6n tv n 07 


838 ± 27 ± 21 ± 17 


77' —> -K^n^jj 27r+27r~27 


824 ± 29 ± 25 ± 18 


77' — > 7r°7r'^77 67 


931 ± 29 ± 21 ± 23 


77' 7r°7r°77 IO7 


837 ± 61 ± 27 ± 17 


Simultaii(H)u.s fit to all // data 


859 i 9 i IS i 20 
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0970597-025 



( a ) 









(b) 



Energy Resolution (%) Angular Resolution (Degrees) 

FIG. 1. Resolution functions of a) energy (in %) and 
b) scattering angle (in degrees) obtained from MC sim- 
ulation in the 7r° analysis. Dashed and solid lines show 
resolution functions measured directly in the calorime- 
ter and achieved using energy-momentum conservation, 
respectively. 
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m(,y~f) (MeV/c^) 
FIG. 3. Fit (solid line) to the 77 invariant mass dis- 
tribution observed in data (points with error bars) in the 
tt" 77 analysis. The signal line shape is obtained from 
the MC simulation, the remaining random background 
is approximated by an exponential. 
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FIG. 2. The efficiency (in %) of the fractional 
tag-energy cut as measured from data. The solid line 
shows a power law approximation chosen to interpolate 
between the efficiency measurements. Events with tags 
scattered at polar angles less than 15° are rejected from 
aU analyses. 
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FIG. 4. Fit (solid line) to the 77 invariant mass dis- 
tribution observed in data (points with error bars) in the 
rj —> 'y'y analysis. The signal line shape is obtained from 
the MC simulation, the remaining random background 
is approximated by the sum of an exponential and a con- 
stant. 
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0970597-026 



^ -in' 




FIG. 5. The distributions for signal tt" events in 
MC (solid line) and data (points with error bars) in the 
Tr" _> analysis. The distributions for events which 
belong to the energy- and track-tagged MC samples are 
shown with dashed and dotted lines, respectively. For 
each interval in data the number of signal events is 
obtained from the fit followed by the background sub- 
traction. The number of MC events is normalized to the 
number of signal tt" events in data. 
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□1(377°) (MeV/c^) 

FIG. 6. Fit (solid line) to the Stt" invariant mass dis- 
tribution observed in data (points with error bars) in the 

^ 67 analysis. The signal line shape is obtained from 
the MC simulation, the remaining random background 
is approximated by a first-order polynomial. 



60 



50 



40 



> 

<i> 

■n 30 h 

en 
c 

<D 

111 20 



10 



850 900 950 1000 

mCjjTrVXMeV/c^) 



1050 



FIG. 7. Fit (solid line) to the Tr^vr^r; 67 invariant 
mass distribution observed in data (points with error 
bars) in the "q' — » 67 analysis. The signal line shape is 
obtained from the MC simulation, the remaining random 
background is approximated by a first-order polynomial. 
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FIG. 8. Fit (solid line) to the Stt'' invariant mass dis- 
tribution observed in data (points with error bars) in the 
77' ^ IO7 analysis. The signal line shape is obtained from 
the MC simulation, the remaining random background 
is approximated by a first-order polynomial. 
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FIG. 9. Fit (solid line) to the Tr^Tr^Tr" invariant mass 
distribution observed in data (points with error bars) 
in the rj —> 7r^7r~27 analysis. The signal line shape is 
obtained from the MC simulation, the remaining random 
background is approximated by a first-order polynomial. 
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m(?j7r*7r") (MeV / c^) 

FIG. 11. Fit (solid line) to the -k^tt^t] 7r+7r~27 in- 
variant mass distribution observed in data (points with 
error bars) in the rj' —> 7r^7r~27 analysis. The signal line 
shape is obtained from the MC simulation, the remain- 
ing random background is approximated by a first-order 
polynomial. 
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FIG. 10. Fit (solid line) to the 7r^7r~7 invariant mass 
distribution observed in data (points with error bars) 
in the rj' — > 7r"'"7r~7 analysis. The signal line shape is 
obtained from the MC simulation, the remaining random 
background is approximated by a first-order polynomial. 
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FIG. 12. Fit (solid line) to the 7r+7r"r7 ^ 27r+27r~27 
invariant mass distribution observed in data (points with 
error bars) in the rj' —^ 2-K^2-K~2'y analysis. The sig- 
nal line shape is obtained from the MC simulation, the 
remaining random background is approximated by a 
first-order polynomial. 
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FIG. 13. Fit (solid line) to the tt+tt'Stt" invariant 
mass distribution observed in data (points with error 
bars) in the 77' — > 7r^7r~67 analysis. The signal line 
shape is obtained from the MC simulation, the remain- 
ing random background is approximated by a first-order 
polynomial. 
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FIG. 15. Distribution of | cos 6**1 in the 77' 
analysis in data (points with error bars) and the MC sim- 
ulation (histogram). The dotted line shows the sin^ 
curve. The prediction of the MC simulation and sin^ 6* 
curve are normalized to the number of data events. 




Signal Photon Energy (MeV) FIG. 16. Measured (points with error bars) and 

FIG. 14. Distribution of signal photon energy in numerically estimated (histogram) differential cross see- 
the T)' p°7 analysis in data (points with error bars) tions for vr" production, 
and the MC simulation (histogram). The prediction of 
the MC simulation is normalized to the number of data 
events. 
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FIG. 17. Measured (points with error bars) and 
numerically estimated (histogram) differential cross sec- 
tions for rj production in the 77 ^ 77 analysis. 
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FIG. 19. Comparison of the results (points) for 
Q^\J-^,^^o (Q^ ) I with the theoretical predictions made by 
Cao et al. ||l^ with the asymptotic wave function (solid 
curve) and the GZ wave function (dashed curve). 



' ' ' I 
□ CELLO 
♦ CLEO 



2.5 




5.0 

(GeV^) 



7.5 



10.0 



FIG. 18. Comparison of the results (points) for 
Q^|J?-"^»^^o (Q^)| with the theoretical predictions made 
by Jakob et al. with the asymptotic wave function 
(solid curve) and the CZ wave function (dashed curve). 
The dotted curve shows the prediction made with the 
CZ wave function when its QCD evolution is taken into 
account. 
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FIG. 20. Comparison of the results (points) for 
Q^\J-ryt.y.,^o{Q^)\ With ths theoretlcal prediction (curve) 
made by Radyushkin et al. M . 
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FIG. 23. 



FIG. 21. The interpolation given by Eqn. |^ (solid 
curve) and the pole- mass parameter fit (dashed curve) to 
our results (closed circles) for \T^,^^o{Q^)\'^ represented Q'^\^-y*wiQ^)\ ^orra (dashed line), 
in the Q'^\T^-.^T,a{Q'^)\ form. 
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Results of the pole-mass parameter fit to 
our results (points) for \T.^*^ri'{Q'^)\'^ represented in the 
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FIG. 22. Results of the pole-mass parameter fit to 
our results (points) for |.7-".y»,y^ (Q^)p represented in the 
Q'^\J-'^'''^-q{Q'^)\ form (dashed line). The solid curve shows 
the interpolation given by Eqn. 
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